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Abstract. The general phenomena associated with sustained resonance are studied 
in this paper in connection with relativistic binary pulsars. We represent such a 
system by two point masses in a Keplerian binary system that evolves via gravitational 
radiation damping as well as an external tidal perturbation. For further simplification, 
we assume that the external tidal perturbation is caused by a normally incident 
circularly polarized monochromatic gravitational wave. In this case, the second-order 
partially averaged equations are studied and a theorem of C. Robinson is employed to 
prove that for certain values of the physical parameters resonance capture followed by 
sustained resonance is possible in the averaged system. We conjecture that sustained 
resonance can occur in the physical system when the perturbing influences nearly 
balance each other. 
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1. Introduction 

Imagine a stable Hamiltonian system tliat is subjected to an external periodic 
perturbation as well as a frictional force that continuously takes energy out of the 
system. If the driving force vanishes on the average, then in general the average behavior 
of the system is characterized by damped motion. Under favorable conditions, however, 
the system could get captured into a prolonged resonance if the driving frequency is 
commensurate with an internal frequency of the system. In this paper, we study the 
dynamics of this phenomenon in the context of a simple model of certain relativistic 
binary systems that are of current interest. 

Relativistic binary pulsars are generally expected to occur in some astrophysical 
environment and the gravitational effect of the surroundings on the binary orbit 
would be tidal in the first order of approximation. The evolution of the binary orbit 
under the combined influence of gravitational radiation reaction and the external tidal 
perturbation can lead to resonance under appropriate conditions [1-6]. Once we average 
over the external periodic tidal perturbation, the binary pulsar is expected to lose energy 
to gravitational radiation damping so that the relative orbit would spiral inward unless 
there is resonance capture. This would come about if the external source deposits energy 
into the orbit over an extended period of time so as to compensate for the continuous 
loss of energy due to radiation reaction. This situation, which persists until the system 
leaves the resonance, can occur if the Keplerian period of the osculating ellipse becomes 
commensurate with the period of the external perturbation. In this way, the binary 
pulsar can be captured into resonance with its gravitational environment. Then, the 
standard data reduction schemes are expected to reveal that during resonance capture 
the semimajor axis of the binary is not monotonically decreasing due to the emission of 
gravitational radiation but is instead oscillating about a fixed resonance value since the 
external perturbation continuously deposits energy into the orbit during resonance so 
as to replenish on average the energy lost to gravitational waves. The system leaves the 
resonance when this delicate energy balance is significantly disrupted. Such resonances 
are well known in the solar system, where the damping effects have a different origin 
(see and the references cited therein). One expects that over time the same kind of 
phenomenon will be observed in a sufficiently large sample of relativistic binary pulsars 
i- 

It should be mentioned that there are transient chaotic phenomena that are 
generally associated with resonance. The transient chaos per se may be difficult to detect 
observationally due to noise and other complications that exist in the data reduction 
process 0. On the other hand, there is numerical evidence for a peculiar chaotic 
effect, namely chaotic transition, that appears to be associated with transient chaos 
and involves a sudden transformation in the rate of orbital collapse That is, the 

system makes a transition from one relative orbit to another that collapses much more 
rapidly. Presumably such a chaotic transition — if it indeed occurs in nature — should 
be detectable by the timing observations of relativistic binary pulsars. Though the 
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unperturbed binary systems in our investigations always involve standard astronomical 
masses on Keplerian orbits, it is possible to extend the analysis to geodesic orbits around 
black holes. Such general relativistic systems involve stable Kepler-like orbits as well as 
unstable orbits. We expect that our analysis of transient chaos could be extended to the 
stable orbits, while the occurrence of chaos involving the unstable orbits perturbed by 
incident gravitational radiation has been demonstrated by Bombelli and Calzetta ||1CI|| . 
These results could be significant for the study of chaotic accretion onto black holes and 
possible chaotic enhancement in the emission of gravitational waves associated with 
active black holes. 

The purpose of the present work is to show that resonance capture can indeed 
occur on average for a simplified model of the physical situation that has been described 
here, since our previous work on sustained resonance has been mainly numerical [2,4,5]. 
Our mathematical model is briefly discussed in the next section. The second-order 
partially averaged system is given in section ^ The geometry of resonance capture for 
the second-order partially averaged system is presented in section |^, and it is shown 
that resonance capture can indeed occur for this partially averaged model under certain 
conditions. The proof of a resonance capture mechanism for the full model equations 
remains an important unsolved problem; instead, we provide numerical evidence for 
sustained resonance in our model in section ^. 



2. Newton— Jacobi equation with radiation damping 

In a series of papers [1-6], we have studied the long-term evolution of a Keplerian 
binary system coupled to a gravitational environment. In fact, the influence of the 
environment has been replaced — for the sake of simplicity and definiteness — by the 
tidal perturbation of an external gravitational wave. More precisely, we have considered 
the simplest gravitational model for a binary system under radiative perturbations. 
In the lowest (quadrupole) approximation, the external wave exchanges energy and 
angular momentum with the binary system but not linear momentum. For the purposes 
of the present investigation, we limit our considerations to a normally incident plane 
monochromatic external wave that is circularly polarized. The amplitude of this external 
perturbation is e, < e ^ 1; it indicates the deviation of the spacetime metric from the 
Minkowski metric as a consequence of the existence of the incident radiation field. In 
practice, e could be extremely small; in fact, direct laboratory detection of gravitational 
waves from astrophysical sources with e ~ 10~^° is the goal of current experimental 
efforts. The relative motion of the binary is planar due to the transversality of the 
external wave and is given by 

— ^ + ^ + = -tKijx' , 1 
at^ P'^ 

where = X\ — determines the relative orbit, the indices i,j G {1,2}, p = \x\, 
and K = Go{Mi + M2). Here the binary system consists of two point masses Mi and 
M2 with the center of mass of the system at rest at the origin of spatial coordinates, i.e. 
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xi{t) = M2 x{t)/{Mi + M2) and a^aW = -Mi x{t)/{Mi + M2). The radiation reaction 
term is given to lowest order by 



7^ 



iGi M1M2 



where p{t) and describe the relative motion in terms of polar coordinates in the 

orbital plane and an overdot indicates differentiation with respect to t. The influence 

of the external radiation on the relative motion is represented in equation (1) by a 

symmetric and traceless tidal matrix /C evaluated at the position of the center of mass, 

^9 r cos fit ± sin fit] 
K. = aQ'^ , . , (3) 

± sm i It — cos i It J 

where fl is the frequency of the external wave and a, which is of order unity, is its 
amplitude. The upper (lower) sign indicates right (left) circularly polarized radiation in 
the transverse-traceless gauge. A general comment is in order here regarding the physical 
realizability of such a perturbation. We can, for instance, imagine another binary system 
with a relative orbit that is circular with Keplerian frequency Q/2. Then, in the rest 
frame of this binary's center of mass and neglecting any dissipation the system emits 
gravitational radiation of frequency Q and the waves that propagate perpendicularly 
to the orbital plane of the circular binary system are circularly polarized; far from the 
circular binary, these waves are nearly planar and of the form under consideration here. 

In the absence of damping, we have shown that if e is sufficiently small, the planar 
relative motion is bounded for all time as a consequence of the Kolmogorov-Arnold- 
Moser (KAM) theorem The physical reason for this confinement is that the 

external wave does not monotonically deposit energy into the binary orbit; in fact, the 
general situation is that energy fiows back and forth between the wave and the binary so 
that on average no net transfer of energy takes place. When the orbital damping is taken 
into account, the general tendency of the binary system is to collapse; however, under 
certain circumstances the system could get captured into resonance. The necessary 
condition for resonance is the commensurability of the Keplerian frequency u and the 
wave frequency Q; that is, relatively prime integers m and n must exist for which 
muj = nQ. Sustained resonance would actually come about if a delicate balance could 
be established between the radiative perturbations. 

It is useful to express the equations of motion in dimensionless form as in our 
previous work [3-5]. To this end, we express all spatial intervals in units of i?o > 
and all temporal intervals in units of Tq > 0, where the scale parameters Rq and Tg 
are connected via hTq/R^ = 1. Specifically, we assume that at some "initial" time the 
lengthscale Rq is the semimajor axis of the binary orbit and 2ttTq is its period. In terms 
of polar coordinates {p,"^), the equations of motion then take the form 

p'^ p'^ P V P 3p / 
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= 25^ i^Vl - 6 ^ + - V sin {2§ t (4) 

Here 6 is the dimensionless parameter that characterizes radiation damping and is given 
by 

' = ^^n^- 

For physical systems, 5 < 2^/^/40 |||]. We find that 5 ~ 10~^^ for the Hulse-Taylor 
relativistic binary pulsar PSR B1913 + 16; this same value for 5 is also approximately 



valid for the relativistic binary pulsar PSR B1534 + 12 |TT . 

To study the general characteristics of system @, it is useful to express this 
system in terms of Delaunay's elements. These are the natural action-angle variables 
for the Kepler problem. At each instant of time t, the state of relative motion {x,x), 
or equivalently (p, 'd, Vp, V^), defines an osculating ellipse that is tangent to the 
perturbed motion at t. The Delaunay variables (£, g, L, G) are closely related to the 
orbital elements of the osculating ellipse. Let a and e be the semimajor axis and the 
eccentricity of the osculating ellipse, respectively; then, P| = a(l — e^), VpV^ = esinv, 
and = 1 + ecos-0, where v is the true anomaly of the osculating ellipse. Moreover, 
let u be the corresponding eccentric anomaly; then, the Delaunay elements are defined 



by 0,[T1 



■u — e sm-u, g 



^-v, L:= a'/^, G:=V^ = L(l - e^^^ . (6) 



Only positive square roots are considered throughout this work. Thus we limit our 
considerations to orbits with positive orbital angular momentum G. 
The Delaunay equations of motion can be written as 

L-^ oL 
9 =e^ + eA7^„ 

G = -e^ + eA7^G, (7) 
dg 

where A = 6/e, TZd, for D G {£, g, L, G}, are radiation damping term^ p, and 
Kxt = ^aa^fi^l^e^ cos {2g t ^t) 

oo ^ 

+ ^-[Kl{e) cos{2gT^t + iye) + K'^{e)cos{2gT^t-ui)]y (8) 

f A misprint in the ^-component of the radiation damping terms in equation (18) of |^ and Hkewise 
equation (A2) of Q should be corrected: In the term proportional to inside the square bracket, the 
factor (1 + 70LV3 - 29G^/2) must be replaced by (70LV3 - 27G^/2). 
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Here 

K^(e) ■.= ^u{A,±B,), 

A, := -4^[2ue{l - e^)Jl{ve) - (2 - e^) J.(z/e)], 

5, := -^(1 - e^fl^ [eJUue) - u{l - e^) J.(z/e)], (9) 

where Jy{x) is the Bessel function of the first kind of order u and Jl{x) = dJ^{x)/dx. 
Finally, we note that A will be viewed as a fixed constant; that is, we will only consider 
perturbations in the (e, 6) parameter space along lines with equations of the form 6 = Ae. 

Let us note that in system (|^ there is only one "fast" angle, namely the mean 
anomaly i, and its frequency is the Kepler frequency {u = L^^) of the osculating ellipse. 
Since Tiext is explicitly time dependent, resonance occurs when muo = nQ, in which case 
L becomes fixed and equal to L^. The dynamics of the system at the (m : n) resonance 
is discussed in the next section. If the system is off resonance, however, the system (J^) 
can be averaged with the result that Tiext averages out to zero and the binary system 
simply collapses due to gravitational radiation damping. 



3. Partial averaging 

The average value of Tiext at the (m : n) resonance is in general nonzero only for n = 1; 
otherwise, the binary system simply collapses on average just as it does off resonance. 
Once the system is captured into a primary (m : 1) resonance, the resonance condition is 
not rigorously maintained. There are in fact deviations that have general characteristics. 
To study these, we let 

i=^t + (j), L = L, + e^l'^V, (10) 

Li 

where and T> are new variables associated with resonance and e^/^ is the corresponding 
small parameter. The average energy exchange due to the external perturbation is 
generally oscillatory 

i'Hcxt) = Tc cos m(f) + Ts siiam(f), (11) 

whereas the damping is unidirectional; in fact, we expect on the average slow oscillatory 
behavior about the resonance manifold. This motion can be described in terms of 
a general damped pendulum-like equation in with torque exerted by the radiation 
damping. 

The equations of motion averaged around the (m : 1) resonance and expressed to 
second order in e^^^ have been given in |^ for the planar Kepler system under radiative 
perturbations and in ^ for the corresponding three-dimensional case. Restricting the 
results of |^ to the circularly-polarized incident wave under consideration here (see 
equations 18-19 of [Q), the second-order partially averaged equations at the (m : 1) 
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resonance are 



■ cos md) + -rrr sm mct 



dL 



dL 



9 
V 



COS m0 + 



dG 



sm mc 



-eV 



— mTc sin m0 + mTs cos m(f) + + "^^^ + Y2^^ 
sinm0 + m^:^ cosm0 + „ ^^^c (,14o + 37e ) 



dL 



dL 



a 



— e 



dg 



cos m(/) + 



sin md) + 



(12) 



(13) 



Here Tc = —/cos 2^, T, = ±/sin2^, and 

f := f{L,G) = -]^amLfL'K^ 

where e = (1 — G^/L^)^/^. The general behavior of the functions K^(e) has been 
discussed in (see figure |I] and appendix A of p). The functions K^{e) and ^^+(6) 
are both negative for < e < 1. In general, the derivative of the function e h-> K^{e) 
is negative for < e < 1 (see p. 107 of ||I|). Moreover, for m > 1, K^{e) is positive for 
< e < Cm, zero for e = Cm, and negative for Cm < e < 1; in fact, 62 — 0.76, 63 ~ 0.85, 
etc. 

It is interesting to note that for incident radiation that is circularly polarized, 
equation ( [TT| ) implies that (Tiext) = — / cos (m0 ± 2^). Introducing a new variable 
6' := mcf) ± 25? and a new slow temporal parameter i := —^t/b, where /i := e^^^ and 
6 := L^/(3m), system ([12|) reduces to the form 

e = V-fi \^V^ - b{mfL ± 2 fa) cos ^ 



± 2/ sin 9 + 



Ab 

G^^" ' 3 
A 



^ 73 37 
V = mbf sin ^ + ^(8 + — + — e^) + /ifeP m/^ sin 6 + 



A 



3L3G5 



;i46 + 376^) 



+ 7e' 



(14) 



where the overdot now refers to the new time t and all relevant quantities are evaluated 
at L = Li,. For example, fi = {df /dL) 



etc. 



It follows from equation ([T^) that the first-order averaged equations reduce to 
G = G^, where is a constant, and T) = 9, where 9(i) is a solution of 

6''-Asin^ = r. (15) 

Here A and r are constants given by 



A = —-amL^K. 



6 
Ab 

Gl 



T3 o 37 A. 

+ y^^ + i2^^)' 



(16) 
(17) 



where e^, 



Gl/LlY^"^. Equation (|l^) describes a mathematical pendulum with 
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Capture Region 




Figure 1. A schematic phase diagram on the phase cyhnder for the simple 
mathematical pendulum with constant torque — A sin 9 — t for the case A > r. 

a constant torque r > 0; the phase portrait of this dynamical system is illustrated 
in figure ^ for A > r. In this case, there is a homoclinic orbit which encloses a 
region corresponding to orbits that are captured into resonance. The region outside the 
homoclinic loop, with the stable and unstable manifolds removed, consists of orbits that 
pass through the resonance. The presence of the perturbation terms in the system (0) 
implies that the pendulum is damped (or antidamped) in the second-order averaged 
equations. That is, averaging the dynamical system to second order results in a 
pendulum equation of the form (|15|) with slowly drifting A and r as well as an additional 
term proportional to ^9. The motion is geometrically described in the next section. 

There are several important phenomena associated with resonance; for instance, 
chaos and resonance capture. Chaotic motion is expected to occur near a resonance 
under a Hamiltonian perturbation. For non-Hamiltonian perturbations, while chaotic 
motions are certainly possible, resonance capture is often the dominant phenomenon. 
At first order, the partially averaged equations at a resonance are generally pendulum- 
like equations with torque. This nonlinear Hamiltonian system can exhibit chaos and 
resonance capture under perturbation. The relationship between the chaos in the 
original dynamical system and possible chaotic effects in the averaged system near 
(m : 1) resonance is rather subtle and beyond the scope of this work. Here, we will 
consider instead the possibility of resonance capture. Indeed, we will give a geometric 
description of resonance capture in system (|T^ and a rigorous proof that resonance 
capture does in fact occur. 
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4. The geometry of resonance capture and a theorem of Robinson 

We start from our model equations in Delaunay variables (|^) and introduce the 
transformation ( [TI1| ) about a primary resonance; then, the resulting system can be 
transformed into the system ( p!^ plus terms of higher order by means of an averaging 
transformation. These higher-order terms may be nonnegligible and dealing with this 
problem is beyond the scope of our work; hence, we consider here resonance capture only 
in the second-order partially averaged equations and not in the full model. System (0) 
is typical of the following class of differential equations 

y = Fo{y, z) + z) + 0{fi^), 

i = /iZi(y,2) + 0(/i2), (18) 

where G and z G M^. In particular, note that z is changing on a slow time scale 
relative to y. For system (|1^, the pair {0,V) plays the role of y, and G, in our case a 
scalar variable, plays the role of z. 

For the geometry of resonance capture, let us make a basic hyperbolicity assumption 
about system (0), namely, let us assume that for each fixed C in some bounded region 
Z of the planar differential equation 

y = Fo{y,C) (19) 

has a hyperbolic saddle point -P(C) with a corresponding homoclinic loop given by a 
solution t y(t, (^). Then, the set 

Mo:={{y,z):y = P{z),zeZ} 

is a normally hyperbolic invariant manifold (consisting entirely of rest points) for the 
unperturbed system (pTSf), i.e. 

y = Foiy,z), z = 0. (20) 

A rough description of this manifold is that nearby orbits are attracted or repelled 
from the vicinity of Mq exponentially fast and the rates of expansion and contraction 
dominate the corresponding rates for orbits on Mq. Of course, in our case where the 
manifold consists of rest points, the second condition is satisfied automatically. The 
normally hyperbolic manifold Mq forms the "spine" of the submanifold M in X 
consisting of the manifold Mq — given as a graph over an open subset of — and all the 
corresponding homoclinic loops as illustrated in figure |^. The manifold M "bounds" 
a region in the phase space for system (pOD, and this region corresponds to the initial 
conditions for orbits that are captured into resonance relative to system (pO]). When 
the system is subjected to perturbation, M is expected to "split", for example, as 
depicted in figure 0. In fact, by a celebrated theorem of N. Fenichel the normally 
hyperbolic manifold Mq persists. If the perturbation parameter /i is sufficiently small, 
then the perturbed system has an "invariant" normally hyperbolic manifold that is 
again given as a graph over a region in M^; in particular, the manifold has stable 
and unstable manifolds. Under appropriate conditions that will be specified below, the 



Sustained resonance 





Sustained resonance 



11 



stable and unstable manifolds of M^, the manifolds resulting from the split homoclinic 
manifold M, form a channel in the phase space where orbits of the phase flow enter the 
vicinity of the resonance (cf. figure This is the geometric description of resonance 
capture for a system such as (|T8|). 

Let us note that because the unperturbed "slow manifold" Mq consists entirely of 
rest points, it is an invariant set for the unperturbed flow. However, the perturbed 
manifold generally has no fixed points; rather, the flow on is characterized by 
a slow drift in the z variable. Thus, is only invariant in the sense that orbits with 
initial conditions on this manifold stay on it until they exit at one of its boundary points. 
Of course, as an effect of the slow drift, an orbit can, after a (perhaps long) sojourn 
near the resonance, enter a new dynamical regime. The sojourn time can be predicted 
from an analysis of the flow on the slow manifold. Thus, our geometric picture gives a 
rather complete description of the phenomena associated with resonance capture. 

To prove that the geometric description of resonance capture obtains for 
system (|18|), we must detect splitting of M and determine the positions of the stable 
and unstable manifolds of M^. This is accomplished by an adaptation of Melnikov's 
method that was first carried out by Robinson [|15|. To state his result, let us choose a 
portion of a plane £ in x that intersects each homoclinic loop in M perpendicular 
to the unperturbed vector field on M. The choice of C gives a smooth choice of initial 
conditions for the unperturbed orbits that correspond to the homoclinic loops. Indeed, 
for each ^, consider the corresponding saddle point -P(C) and the unperturbed solution 
t yo(t,C) such that yo{0,C) corresponds to the intersection point of £ with the 
homoclinic loop at -P(C)- Let d{^,C) denote the distance between the intersections of 
the stable and unstable manifolds of as measured along C in the intersecting plane 
parallel to the y coordinate plane and passing through the point (P(C),C)- Also, for 
two vectors a,b ^ M^, let a A 6 := 0162 — ct2^i- The following theorem is proved in |15]: 

Theorem 4.1. If system ( fl^ j satisfies the basic hyperbolicity assumption and the vector 
field given by y -fo(?/)C) ^■^ divergence-free for each fixed ^ E Z, then the splitting 
distance is given by d{fi, C) = At'^i(C) + 0(/^^)) where the leading-order coefficient is given 
by the integral 

all functions in the integrand are evaluated at {yo{t, C); C); <^iT'd the function 
dz 

t^gj^{yo{t,C)X) 

satisfies the variational initial value problem 

r = Zi(yo(t,C),C), Y{o) = o. 



The following immediate corollary is also stated in [ p!5|] : 

Corollary 4.2. If di{(^) > for all C, E Z, then there is resonance capture. Moreover, 
the only way a captured solution can leave the vicinity of the resonance is to have its 
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second component, z, leave the region where di{(^) > 0, a process that occurs on a 
timescale of order 1/yU. 

The existence of resonance capture means that a capture channel opens up as in figure ^ 
and therefore there is sustained resonance. 

Let us now formulate and prove a version of another result of Robinson that we 
use to apply theorem ^?T| to the pendulum-type systems that arise by partial averaging 
at resonance. Recall that system ( p!4|) is derived from the averaged physical model in 
part by a reversal of time. Thus, it is convenient to have a result that implies "release 
from resonance" for system (|14|) so that it will imply capture into resonance for the 
original averaged model. Of course, the phrase "release from resonance" means that 
the stable and unstable manifolds of the normally hyperbolic manifold Mq split in the 
direction that is opposite to that required for capture; in effect, orbits leave the vicinity 
of the resonance in a channel bounded by the unstable manifold. 

Proposition 4.3. Consider the system 

if = V + ij,p{(p,v,w), 

V = f{w)sm(p + g{w) + ij,q{ip,v,w), 

w = fi r{{p,v,w), (22) 

where {ip, VjW) E M. x M. x M^, yU > 0, and all indicated functions are smooth. If there 
is a point w G such that 

f{w)>0, g{w)>0, f^<l, (23) 

then the unperturbed system 

= f , V = f{w) sirup + g{w) 

has a homoclinic loop that encloses a region R{w) in the {(p,v)-plane. Suppose that, in 
addition, t ^— {(f{t,w),v{t,w)) is a solution on the homoclinic loop with 

f'{w)r{ip{t,w),v{t,w),w)<0, (24) 

r{ipit,w),vit,w),w) > 0, (25) 

where a prime denotes differentiation with respect to w, and for {^p^v) G R{w), 

p^{(p,v,w) + qy{(p,v,w) >0. (26) 

// at least one of the inequalities (j^j, (IH^j or (j^^ is strict, then Robinson's integral 
= d,i in is negative atw. In particular, the perturbed configuration of stable and 
unstable manifolds corresponds to "release from resonance". 

Proof. The first statement of the proposition, about the existence of homoclinic loops, 
is easily proved by elementary phase-plane analysis. Let s denote a new time-like variable 
given by 

ds = {f{w))'/'dt, 
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such that s = at t = and let V := f /(/(u;))-^/^. In these new variables, system 
has the form 

ds (/M)V2 ' 

dV . g{w) f q{^,{f{w)Y/'V,w) 

— = sm + + /i — 

ds f{w) V f{w) 

-^ifiw))-lf'{w)Vr{v, (/M)^/V, 

dw _ r{ip,{f{w)y/^V,w) 



w 



ds " ifMy/' 

Robinson's integral Jr is then given by 



Ik-- 
where w 





g/f 



w, 



A 



V 

simp + g/f 



ds, 



^ .— dw/dfi. We find that Jr = Ji + /2 + /s, where 

/oo 
[pf-'/' ismif + g/f)-Vq/f] ds, 
-oo 

Vig/f)'w,ds. 



oo 
oo 



Also, let us note that w is constant on the separatrix loop and let us take the solution 
of system (p2D such that the partial derivative vanishes at s = as required in the 
statement of theorem ^T^. The first integral Ji can be rewritten as follows: 

/oo 
sin cp + g/f) - vf-'/'q/f f'/' dt 
-oo 

= 7/ [p{fsmip + g)-vq]dt. 

J J — oo 

Under the hypothesis (|23|), the homoclinic loop F in the {ip, f )-plane is parametrized with 
negative orientation relative to the usual orientation of the plane. Also, let us denote 
by R the bounded region with boundary F in the {ip, f )-plane. By Green's theorem and 
the inequahty (p6D, we have that 
1 



^-7 



/ 



(pi) — qtp) dt = —— I pdv — qdip = —— I {p^ + q^) d(pdv <0 



f 



R 



It follows from hypothesis (|2^) that /2 < 0. For the integral Is, let us observe that if 
we take the parametrization of F such that V{0) = 0, then sV{s) > 0; that is, V{s) 
has the same sign as the variable s in the integration. Also, because w is constant for 
the integration, ju = for the functions in the integrand, and hypothesis (^) holds, we 
have the inequality 



d_ 

d^ 



w, 



gydw^ 



fJ ds 



^l=o 



(f)'rr"^>o. 
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Thus, using the fact that vanishes at s = 0, note that {g/f)'w^ has the same sign 
as the variable s for the integration. We conclude that 

over the range of the integration. Thus, I3 < 0. 

Finally, note that if inequality (|26|) is strict, then Ji < 0. Likewise, if inequality (0) 
(respectively inequality (|25|) ) is strict, then I2 < (respectively I3 < 0). Hence, 
-^R = -^1 + -^2 + -^3 < and if at least one of the inequalities (p^, (p5D, or (p6D is 
strict, then Jr < 0. □ 



To apply proposition let us consider the second-order averaged equations for 
the normally incident left circularly polarized gravitational wave expressed in decoupled 
form; namely, system (0) with the "lower" signs. Also, let us note that this system is 
a special case of the abstract system in proposition |4.3| once we set 6 = (p, V = v, and 
G = w. Recall that if > G > 0, then < e < 1 and K^{e) < 0. Hence, we have the 
desired inequality f{L^, G) > once we choose a > 0. In fact, we set a = 1. We claim 
that the hypotheses of proposition [4.3| are satisfied under the following two conditions: 

1. There is a subinterval Qi of the interval (0, L^) such that 

73 2 37 4 , ^ 
-e + — < 1 



m/G7 V 3 12 
whenever G E Qi. 

There is an interval Q2 ^ Gi such that the derivative of the function 

1 /„ 73 2 37 4 



(l_e2)7/2jV 3 12 
is negative on the image of the interval Q2 under the transformation 
e = (1 - G^Liy/^ 

Condition 1 ensures that hypothesis (p3D is satisfied for each G E Qi. Recall that 



dedG Ge^ ' de' 
and since dK"^{e)/de < for < e < 1 we have that df/de > 0. Thus, /g < for 
G E Qi. Turning to the function r in the statement of proposition [4.3| , note that 

Ab 

r:=-26/sin0 + -^(8 + 7e2). 



By the first implication of proposition [4.3| , it follows that if G G ^1, then the system 
e = V, V = mbf sine + [8 + + ^e'^ 

has a homoclinic orbit where the corresponding saddle point has coordinates {6, T>) = 
{Oi,, 0) for some 6'^ in the interval (37r/2, 2tt) as long as 6' G [0, 27r). At the boundary of 
the interval Qi, where 

A f 73 „ 37 4 



mfG^ V 3 12 
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the corresponding saddle point with coordinates {9^, 0) and center with coordinates 
(Svr — 9^,0) coalesce at (37r/2,0); that is, these rest points disappear in a saddle node 
bifurcation. Thus, there is a subinterval Q2 ^ Gi such that if G E Q2, then the region 
R{G), enclosed by the homoclinic orbit, contains only points whose first coordinates 9 
are in the interval n < 9 < 27r. In particular, if G G Q2, then sin 6' < and therefore 
r > 0. Hence, rfc < and hypothesis (^) is satisfied. 

To verify hypothesis (^), we must show that if G E Q2, then the derivative of the 
function 

A 73 2 37 4 

G^ — 8 + — 6^ + — 
m/G7 V 3 12 

is positive; however, this requirement follows immediately from condition 2. For 
hypothesis (^), let us determine the sign of the divergence of the vector field V given 
by 

{9,V)^ ^-(^^ + b{2fG-mh)cos9y bV (mh sin 9 + ^^{U6 + 37 e'' 



Indeed, using the inequalities sin 6' < and fc < 0, it follows that if G G Q2, then the 
divergence of V is the positive quantity 



6(2/Gsin^ + ^^(146 + 37e2; 



This completes the proof of our claim. 

To show the desired phenomenon of release from resonance, we will verify 
conditions 1 and 2. In order to determine the interval Q2, let us note that our 
(unperturbed) pendulum-type planar system has the form 

9 = V, V = Xsin9 + T, (27) 

where A > 0, r > 0, and r/A < 1. System (27) is Hamiltonian with total energy 

H = -V^ + Xcos9 -t9. 
2 

Denote the coordinates of its hyperbolic saddle point by (^^^,0), and recall that 
37r/2 < 9^ < 271. Also, note that the homoclinic loop is given by the graphs of the 
functions 

V = ±2^/\X cos 9^ - t9^ - {X cos 9 - T9) f^ . 

Proposition 4.4. For system ^^T), there is a point 9 such that 9^, — 2n < 9 < 9^^ and 

X cos 9 — t9 = X cos 9^, — t9^,. 

The point {9, 0) is not a rest point. Moreover, the equation 

X cos 9 — t9 = X cos 9^, — t9^, 

has no solution for 9 > 9^,. That is, the homoclinic loop at (6'^,0) crosses the 9-axis at 
9 = 9 and encloses that portion of the 9-axis given hy 9 < 9 < 9^,. Also, 9 > tt if and 
only if k := t/X is such that 

l + (l_fc2)l/2 _ 

^ ; h Sm k < TT. 

k 
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Proof. Define the function h by 

h{e) = Xcose -t9-{X cos - rO^) 

and note that it has relative maxima at 6 = 6^, and 6 = 6^, — 27t. Also, h{6^) = and 
h{6^ — 27t) = 2'KT > 0. It follows that there is a point 6 in the interval indicated in 
the statement of the proposition such that h{6) = 0. Moreover, this point is not a rest 
point. Indeed, the only rest point in the open interval (0^ — 27r, 6^,) is the point 371/2 — 6^,, 
and it corresponds to a negative relative minimum of h. If ^ > 6'^, and h{6) = 0, then 
there must be some relative maximum of h that is nonnegative. However, all the relative 
maxima that exceed 6^ have the form 6^ + 2'n'N, where is a positive integer, and h 
takes negative values at all such points; therefore, we reach a contradiction. Clearly, 
6 > Tc if and only if /i(vr) > or 

t{9^ - tt) - A(l + cos 9^) > 0. (28) 

Recall that sin^^ = — r/A and cos6^ > 0. Thus, with k := r/A, inequality (|28|) is 
equivalent to the inequality 

1 + (1 - P)V2 

^ ; h Sm ^ < TT, 

k 

as required. □ 

In summary, let us recall that /(e) = mbf{Li,,G), where e = (1 — G^/LlY^^ and 
f{Li,,G) = —{m/2Ll)K^{e) by equation (p!3|). For incident left circularly polarized 
radiation, we have that /(e) > and /'(e) > 0; in fact, 

I^-m/ X _ (5m + 2)m'" ^ rn+2 , r\(prn+A^ 

2™+i(^ + 2)! ^ ^' 

and for < e < 0.5 the first term in this Taylor expansion is accurate to about 10%. 
Then, all the hypotheses of proposition are satisfied provided that there is some e 
such that < e < 1 and 

/ N , / N A 1 / 73 r, 37 .\ 

« - 7^(13^ + +12^) 

(ii) k'{e)<0, 

H ^(^^ +sm-^(fc(e)) <7r. 

One can show by inspection that condition (iii) is satisfied provided that fco < k{e) < 1, 
where k^ ^ 0.725. We therefore choose k{e) = 0.9 in agreement with condition (i). 
Moreover, one can directly verify that k' (e) < if we assume e = 0.5 and m = 2. It is 
interesting to note that the eccentricities of the relativistic binary pulsars PSR B1913+16 
and PSR B1534+12 are ^ 0.6 and ^ 0.3, respectively. We choose U = (4/3)^/^ so 
that Q ~ 1.3 at the (2 : 1) resonance that we wish to consider. Finally, A is uniquely 
determined from k{e) = 0.9 once we recall that a = 1 and ^^^(0.5) ~ —0.009; hence, 
A ^ 5 X 10"^. 
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5. Discussion 

We have shown that a set of parameters exists such that for incident left circularly 
polarized radiation the averaged equations permit sustained resonance. The essential 
requirement for sustained resonance in our approach is that the Robinson integral be 
nonzero. Our objective has been to prove sustained resonance in the second-order 
averaged equations for our model. However, we expect that Robinson's integral is 



generally nonzero without the restrictive assumptions of proposition [4.3| . Indeed, it 



should be clear from the proof of proposition that Jr can be negative without the 
requirement that each of the three integrals Ji, I2, and Is be negative. For this reason, we 
expect that resonance capture would generally be possible for either circular polarization 
state of the incident radiation. Moreover, circular polarization is not essential for 
sustained resonance as demonstrated by our previous numerical work on resonance 
capture that did not involve circular polarization [2,4,5]. Indeed, the mathematical 
theory of sustained resonance developed in this paper owes much to the simplification — 
i.e. the reduction from four to three phase-space dimensions in the system (|14D — brought 
about by the assumption that the normally incident radiation has definite helicity. 

The proof of the existence of sustained resonance involves the second-order partially 
averaged system. We conjecture that the result can be extended to the original system. 
To provide evidence for this conjecture, it is interesting to use the set of parameters 
obtained in the previous section as part of the initial conditions needed to integrate the 
original system numerically. To this end, we integrate system (^) with initial conditions 
(p, ^9, Vp, V^) = (1, 0, 0.5, l),n = 2(3/4)3/2, e = IQ-^, and 5 = eA = 5 x 10"^. Thus the 
initial orbit has L = (4/3)1/2, e = 0.5, G = 1, and ^ = 57r/3- 3^/^/2, so that n< 9 < 27i 
as required. The periastron for this initial osculating orbit, which has true anomaly 
V = it/2, occurs at ^ = — 7r/2; hence g = — 7r/2. In fact, this initial orbit has been 
used frequently in our previous work [3-6]. We find that with these initial conditions 
the system is indeed in sustained (2 : 1) resonance for either circular polarization state 
as demonstrated in figures 4 and 5. Figure | can be produced as follows: System (^) 
is integrated backward from the initial point for normally incident radiation that is left 
circularly polarized until the system leaves the resonance and a final point of integration 
is stored; then, the graph is made by integrating forward in time from this point. It 
is interesting to note that the binary orbit collapses rather rapidly following exit from 
resonance when the orbit is highly eccentric. This is explained by the fact that the rate 
of emission of gravitational radiation by the binary orbit is proportional to (1 — e^)"''/^; 
hence, as e ^ 1 the loss of orbital energy to radiation leads to orbital collapse ||^. 
Figure |^ can be produced in a similar way and illustrates the same situation for right 
circular polarization. The possibility of resonance in this case can be seen from the 
first-order averaged equations (|T5|)-(P!7|); that is, ^^^(0.5) ~ 0.923 and hence |A| > r, so 
that the situation depicted in figure 1 holds except that the direction of some arrows 
should be reversed. It is clear from figures 4 and 5 that evolutionary dynamics of the 
orbit while trapped in resonance depends sensitively upon the helicity of the incident 
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Figure 4. Sustained resonance for the relative orbit in a binary system with initial 
conditions {p,'&,'Pp,V^) = (1,0,0.5, 1) in (2:1) resonance with a normally incident left 
circularly polarized gravitational wave. In the graphs, produced by backward and 
forward integration as explained in the text, this initial point corresponds to f « 9500. 
Here a = 1, e = 10'^, 6 = 5x 10"^, and Q = 2(3/4)3/2 « 1.299. The top, middle, 
and bottom panels depict, respectively, L, G, and e versus time. The rate of orbital 
decay after the resonance is generally different from that before the resonance, since 
during resonance capture the orbit exchanges angular momentum with the incident 
radiation such that the orbit has a different eccentricity when it exits the resonance. 
The duration of the resonance is w 15/e time units, in agreement with corollary 4.2. 
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Figure 5. Same as in figure ^except that the incident radiation has positive hehcity. 
The duration of the resonance is ~ 34/e time units; this is consistent with coroUary 4.2. 
The response of the system is more pronounced in this case, since the sense of the 
orbital angular momentum is the same as the helicity of the incident radiation. 
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radiation. For example, the orbit loses angular momentum while in resonance and 
becomes highly eccentric when the incident radiation has negative helicity as in figure 4, 
while the opposite situation holds in figure 5 for positive helicity incident radiation; this 
circumstance may be naturally interpreted in terms of the (negative or positive) angular 
momentum carried by the (negative or positive helicity) radiation and deposited into the 
binary orbit. The fact that the amplitude of oscillations around resonance is generally 
larger when the sense of the helicity of the radiation is the same as the rotation of the 
orbit is in agreement with our previous work on the Hill system. In contrast to the Hill 
system, the present case involves dissipation; nevertheless, the analysis of the amplitude 
of oscillations around the resonance is very similar to that carried out in section 6 of 
due to the fact that r < |A|. 

Figures 4 and 5 also illustrate the complex structure that is usually associated 
with a higher-order resonance [^; furthermore, the presence of chaos is strongly 
indicated in our numerical work. Chaos is especially noteworthy in figure 5; therefore, 
the qualitative significance of our numerical results should be emphasized. We have 
performed numerical experiments with e and 6 twice as large as those in figures 4 and 5. 
In this case, we find in the experiment corresponding to figure 5 that the (2:1) resonance 
appears to have significant overlap with the (3:1) resonance. Let us note that we expect 
e to be many orders of magnitude smaller in physical situations of interest; however, we 
have taken e to be 10~^ in figures 4 and 5 for the sake of simplicity. 

Do all relativistic binary orbits monotonically decay by gravitational radiation 
damping? There is evidence for such energy loss from the timing analysis of only 



two pulsar systems |TT|]. However, such a system could in principle fall into sustained 
resonance with an external source; during resonance lock, the orbit would not decay 
on the average. The orbital decay would resume once the system is released from 
resonance. It remains to see — as more binary pulsars are discovered and their timing 
data accumulate — whether sustained resonance in fact occurs in relativistic binary 
systems. 
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